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ABSTRACT 

We present observations of the spiral galaxy NGC 7331 using the Sub-millimetre 
Common User Bolometer Array (SCUBA) on the James Clark Maxwell Telescope. 
We have detected a dust ring of 45 arcsec radius (3.3 kpc) at wavelengths of 450 and 
850-^tm. The dust ring is in good correspondence with other observations of the ring in 
the mid-infrared (MIR), CO and radio-continuum, suggesting that the observed dust 
is associated with the molecular gas and star formation. A B-K colour map shows an 
analogous ring structure with an asymmetry about the major axis, consistent with 
the extinction being produced by a dust ring. The derived temperature of the dust 
lies between 16 and 31 K and the gas-to-dust ratio between 150 and 570, depending 
on the assumed dust emission efficiency index (/? =1.5 or 2.). 
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1 INTRODUCTION 

Modelling of COB-E-DIRBE observation of our Galaxy has 
shown that dust at T <22 K associated with neutral and 
molecular gas is responsible for most of the emission (55 



85 per cent) in the far-infrared (FI R) ( Sodroski et al. 1994 
Sodroski et al. 1997] ) . Reach et al. ( |l995| ), analysing COBE- 



FIRAS spectral data, find that Galactic emission between 
100- /^m and 300- /^m is predominantly from dust at T = 16-23 
K. This range of temperatures is in agreement with theoreti- 
cal calculations for clas sical grains heated by the interste llar 
radiation field (ISRF) ( |Disney, Davies fc Phillipps 198S| ). 

In external galaxies, dust temperatures have been 
mainly measured using IRAS fluxes at 60-/im an d 100-/im 



is because of the strong dependence of emission on temper- 
ature (oc T A +' 3 ,I3 = 1 - 2). 

Estimates of the temperature of the main dust compo- 
nent depend on th e availa bility of data at A >100-/^m. Re- 
cently, Alton et al. (1998a) have measured a (T) ~ 20 K for a 
sample of 7 spiral galaxies, using /SOPHOT 200-/im maps 
combined with HiRes IRAS maps at 100-/im, 10K lower 
than the temperatures measured using IRAS HiRes 60-/im 
data. The mean gas-to-dust ratio is then ~ 230, closer t o 
the canonical Galactic value of 160 (sodroski et al. 1994). 



Newly available FIR maps, despite their low resolution (117 
arcsec for KOPHOT data) have revealed that the spatial 
distribution of dust is more extended than that of the IRAS 



IRAS can only trace warm dust at T 30 K (Devereu: 



and K-b and (stellar) emission (Alton et al. 1998a; 



fc Young 1QQ0 ) This is mainly associated with star fnrm- 

ing regions, as illustrated by t he strong correl ation between 
Ha and IRAS-FIR emission ( frPevereux 1995 ). Gas-to-dust 



al. 1998) 



Davies et 



ratios derived f rom IRAS are a factor o f 10 larger than the 
Galactic value ( Devereux & Young 199C ) . Assuming that the 



Galaxy is typical of spirals, the discrepancy can be explained 
by a colder dust component (T ~ 15 K) that contributes to 
the bulk of the mass, with out adding a lot of flux in the 
window observed by IRAS (Devereux & Young 1990). This 
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In this letter we describe sub-millimetre (submm) ob- 
servations of NGC 7331 with SCUBA. Maps produced by 
SCUBA have a far higher resolution than J50PHOT and 
HiRes IRAS long wavelength observations and give us de- 
tailed spatial information about the extent of cold dust. 
NGC 7331 is an Sb g alaxy (| dc Vauc ouleurs et al. 1991 ) at a 
distance of 15.1 Mpc (Hughes et al.1998), giving an apparent 
scale of 73 pc arcsec - . 

NGC 7331 shows a ring structure in CO emis- 
sion ( Tosaki & Shioya 1997; von Linden_ej i _aL_J : 996[jYojnig 
fc ScoviTlc*T?)82| ), in radio-continuum (Cowan, Romanishin 
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JCMT 

of the Particle Physics and Astronomy Research Council of the 
United Kingdom, the Netherlands Organisation for Scientific Re- 
search, and the National Research Council of Canada. 



fc Branch 19941) and in MIR images from ISO CAM (6.75 



^m< A <15-/wn) ( (Smith 1998; ). Smith (J1998J) also shows that 
these observations are in good correspondence with each 
other, indicating the presence of a molecular and dust ring of 
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radius w 45 arcsec (3kpc 



Ther e is also a good match with 



Hq+[NII] features (P oggel989h. The central region (radius 



2 arcm in) is partially depleted of neutral hydrogen ( Bosnia 



1981) whilst the molecular gas contributes to 7 per cent of 
the tot al ga s mass inside the optical radius (Devereux & 
Young L990). A major axis profile at 100-/im obtained with 



the Kuiper Airborne Observatory (KAO) (Smith & Harvey 



1996j) sfiows a flat-topped FIR distribution (FWHM » 35 
arcsec), but very little structure can be seen in the HiRes- 
IRAS maps or in the 200-/im TSOPHOT observations of 
Alton et al. (|l998a| ) [FWHM ^ 50 arcsec). 

The ring may be the result of long-term dynamical evo- 
lution driven by a bar instabilit y, as p roposed by von Linden 
et al. ( [l996|) . Tosaki & Shioya ( 1997|) relate the ring to the 
post starburst status of the galaxy: the central gas is con- 
sumed in the long evolution after the starburst event. 



2 OBSERVATION AND DATA REDUCTION 

NGC 7331 was observed at JCMT with SCUBA on 1997 Oc- 
tober 20, 22, 24. SCUBA consists of two bolometer arrays 
able to simultaneously image a region of sky of about 2.3 ar- 
cmin in diameter: the short-wavelength array, optimised for 
observing at 450-/im (91 elements) and the long-wavelength 
array optimised at 850- /mi (37 elements). To fully sample 
the field of view with both the arrays, the secondary mirror 
is moved in a 64-point jiggle pattern. For each jiggle position 
the integration time is 1 s. Simultaneously, the secondary 
mirror is also chopping with a frequency of 7 Hz to remove 
the sky background. Every 16 steps of the jiggle pattern, the 
telescope nods, to remove slowly varying atmospheric gradi- 
ents. We used a chop throw of 180 arcsec perpendicular to 
the major axis of NGC 7331. 

We frequently determined the transparency of the at- 
mosphere during each night by measuring sky emission at 
several elevations. Sky conditions were stable during most 
of the observing run with zenith optical depths of T450 = 
0.6 — 0.7 and rgso = 0.13 — 0.15. During the first night at- 
mospheric opacity was higher, with T450 > 3 and rgso =0.5. 
The telescope pointing was checked every hour against a 
bright point source close to our target: rms pointing errors 
were « 3 arcsec in both azimuth and elevation. Several im- 
ages of Uranus were taken for photometric calibration. 

Data reduct ion was carried out using the STARLINK 
package SURF fljenness fc Lightfoot ll)^ |Sandell 1997| ). 
After subtracting the off-source signal, images were flat- 
fielded and corrected for atmospheric extinction. Noisy 
bolometers were masked. Each image was corrected for sys- 
tematic noise variations using stable bolometers that ap- 
peared free of source emission. Spikes from transient detec- 
tions were removed by applying a 3<j clip. 

A calibration constant for each night was computed 
from the Uranus maps. Comparing data for each night we 
derived a relative error in calibration of 14 per cent and 7 
per cent, for 450-/im and 850-/im respectively. Uranus im- 
ages were also used to measure the beam size and the con- 
tribution of side lobes. The measured HPBW is 10 arcsec at 
450-/im and 15 arcsec for 850-/im. The side lobes can intro- 
duce quite large systematic errors: for example, 20 per cent 
of the emission of a point source goes in the side lobes for the 
long wavelength and 60 per cent for the short wavelength. 



In the following, values of integrated flux are corrected for 
the side-lobe pickup. 

A total of 15 images at 450-/im and 19 at 850-/im were 
re-sampled to an equatorial reference frame, using pixels of 
3 and 5 arcsec, respectively. The final map covers an area 
of 3.3 arcmin x 5.8 arcmin. The total on-source integration 
time for the central part of the image is 5700s and 8300s. 
Maps are shown in Fig. contours with S/N > 3 are over- 
plotted. The short wavelength image clearly shows a poorer 
signal to noise, mainly due to the higher emissivity of the 
sky at 450- /im compared to 850- /mi: the sky noise is 80 mjy 
beam -1 at 450-/im and 8 mjy beam -1 at 850-/im. 

In this work we have also used B and K-band observa- 
tions. Images were obtained from the Skinakas observatory 
Crete (B) and at the Wyoming Infrared Observatory (K). 
The optical an d NIR data were reduced and calibrated in the 
standard way ( [Trewhella 1997| ). The PSF has a FWHM=1.5 
arcsec for the B-band image and 2.5 arcsec for the K-band. 
Optical images were aligned in a RA/Dec frame using field 
stars in the HST guide star catalogue. The peak emission in 
B-band image occurs at RA=22 h 37 m 04 s , Dec=34° 24' 55 s 5 
(J2000); these coordinates are within 1.5 arcsec of th e opt i- 
cal and radio emission centres given by Cowan et al. ( 1994 ) . 
SCUBA images were aligned by assuming the B peak as the 
centre of the frame. Images in B-band and B-K colour are 
shown in Fig. |l| The B-K colour image has a photometric 
accuracy of 0.2 mag. 



3 RESULTS & DISCUSSION 

Our 850-/im image (Fig. [lj top-right) clearly shows a ring 
structure, of w 90 arcsec x 30 arcsec (radius~ 45 arc- 
sec, corresponding to 3.3 kpc). This submm ring matches 
MIR observations, and therefo re CO , radio-continuum and 
Hq+[NII], as shown in Smith jl998[ ). 

The brightest parts of the ring are on the north and 
south sides, where intensity peaks at almost the same value 
(95 mjy beam -1 , S/N fa 12). The east and west side of the 
ring are less br ight a nd appear clumpy. As suggested by von 
Linden et al. (1996) for the CO ring, this morphology may 
well be due to the inclined view. There are also two smaller 
structures attached to the north and south sources (seen in 
our image at positions [5,60] and [20,-50], respectively). The 
south one, the brightest, is associated with the origin of one 
of the spiral arms in our optical image (Fig. bottom-left). 
There is no e vidence for a central source, contrary to the 
MIR images flSmith 1998; ). The ring is also visible in the 
450-/im image (Fig. |l|, top-left) despite the poorer signal. 
When smoothed to the long wavelength resolution, the same 
structures can be seen. 

Radial surface-brightness profiles were produced for the 
850- /im and the smoothed 450-/im images, averaging over 
elliptic annuli. We have adopted a disk inclination o f 74° 
and PA =167° (iGarcia- Gomez & Athanassoula 199l|). The 



two profiles can be followed out to 150 arcsec from the cen- 
tre (both images with S/N > 1), and they appear similar 
(Fig.§. 

Using this last aperture, we have computed the inte- 
grated fluxes: we obtain a flux of 13.4 Jy for 450-/im and 
1.9 Jy for 850-/im. The ratio between the fluxes is therefore 
7.0 ± 1.0, where the error is mostly due to the uncertainties 
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Figure 1. 450-/im (top-left), 850-/^m (top-right) B-band (bottom-left) and B-K (bottom-right) maps of NGC7331. For submm images, 
contours are overlayed at integer values of the S/N ratio, starting from S/N=3 (for the 450- /xm image this corresponds to isophotes every 
80 mjy beam" 1 starting from 240 mjy beam -1 ; for 850-/im, isophotes every 8 mjy beam -1 starting from 24 mjy beam -1 ). Contours 
at 850-^tm are also overlayed on the B-K image. 



in the calibration. From this ratio the average dust temper- 



ature can be computed, if the dependence of dust emissivity 
Q em on wavelength A is known. It is generally assumed that 
Q em oc A -,s , with (3 increasin g from 1 to 2 go ing from the 
MIR to the FIR and submm ( |Hildcbrand 1983|). F its of the 
Galactic FIR-submm spectrum (Reach et al. 1995) suggest 



a value f3 = 2 for A >200-/im. Millimetr e/submm obser - 
vations of the Galactic plane and cirrus (Masi et al.1995) 
emissions are better described by /3 = 1.5. Using (3 = 2 we 
derive T =16±3 K (T =31±10 K for (3 = 1.5). The tem- 
perature obtained using f3 — 2 agrees very well with values 



for high latitude dust emission in the Galaxy (Reach et al 



1995). It is interesting to note that a grey-body model with 
(3 = 1 cannot be fitted to the data. This is also confirmed 



by SCUBA observations of NGC891 (Alton et al. 1998b) 



The dust temperature was also derived using th e KAO 
100-/im observations presented by Smith & Harvey (1996). 



The image at 850- /im was smoothed to the KAO resolution 
(31 arcsec x 41 arcsec) and a profile along the major axis 
was produced, with the same orientation as the one shown 
in their Fig. 5. In the central 100 arcsec, where the KAO 
profile appears flat, the average flux density per bolometer 
at 100pm is 29 Jy. In the same region the average value 
for the 850-pm flux is 60 mjy beam -1 . From the ratio we 
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Figure 2. Elliptically averaged surface brightness profiles for 850- 
/im and 450-fim. Isophote geometrical parameters are described 
in the text. Error bars rapresent the standard deviation of the 
mean inside each elliptical annulus. 
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Figure 3. Correlation between pixels in a B-K smoothed image 
and 850-/im image, inside a 100 arcsec x 50 arcsec region centred 
on the galaxy. Filled circles refer to the far side of the galaxy 
(i.e. east side) and open circles to the near side (i.e. west side). 
The ring structure evident in 850-fim image (Fig. [lj) correspond 
to fluxes bigger than 60 mJy beam -1 . Each data point has a 
random error of «0.05 mag. 



derived T =23K±2 for j3 = 2 (T = 29K±3 for f3 = 1.5). 
Alton et al. ( |l998a| ) fluxes at 100-^m and ISO 200-Atm give 
values T =17K for j3 = 2 and T =19K for f3 = 1.5. 

The correspondence between the submm ring and other 
observations suggests that the dust emission is associated 
with the molecular gas. Indeed, the values of dust tempera- 
ture we have presented here are consistent with those of the 
dust co mponent associated w ith the molecular gas in the 
Galaxy (sodroski et al. 1997). 

The B-K map is shown in Fig. [I], bottom-right, with 
the 850-/im contour superimposed. A ring structure can be 
seen in the colour image, in excellent correspondence with 
the submm image. The ring is delimited on the east and the 
west side by regions with a redder colour with respect to the 
central one. The colour on the west side is also redder than 
on the east. The asymmetry and the central gap, fi rst ob - 
served in NIR colors by Telesco, Gatley & Steward (1982), 
can be explained with a bulge obscured by an inclined dust 
ring: on the west side, the closer one, most of the bulge 
is behind the dust ring so that bulge light is more extin- 
guished and therefore redder than on the east side, the far 
one. Because of the relative lack of dust in the middle of the 
structure, bulge light from the central region appears bluer. 

To compare B-K with 850-/im emission, the colour im- 
age was smoothed to a resolution of 15 arcsec and rebinned 
on 5 arcsec pixels. We have then plotted B-K vs 850-^tm for 
each pixel in an elliptic region of ~ 100 arcsec x 50 arcsec 
(Fig A clear trend is evident: for the same 850- /im flux, 
the closer part of the galaxy is systematically redder than 
the far one, by 0.5 mag. This is true both for the regions 
with high flux (i.e. high dust column density) that define 
the ring and for regions external to the ring. Bulge light ex- 
t inguis hed by a more extended dust disk may be the reason 



of a possible misalignment of the colour and submm images, 
which we estimate to be at most 2.5 arcsec. 

B-K is m ore sensitive than optical colours to the pres- 
ence of dust (Block et al. 1994), because of the high ratio 
between extinct ion in B and K -band (« 14, using a Galactic 
extinction law ( Whittet 1992j ) ) . Even so, deriving the dust 
content from the colour image is quite difficult: first an a pri- 
ori knowledge of the intrinsic colour distribution is required, 
to quantify the reddening; second, the relative distribution 
of dust and stars must be known, to infer, by means of a ra- 
diative transfer model, the relation between the reddening 
and the optical depth (or the dust column density). These 
two steps are obviously coupled together (See for example 
Xilouris et al. 1997, 1998). 

We can easily obtain a lower limit to the optical depth 
of the ring. We assume the thickness of the dust structure 
to be negligible so that dust acts like an internal screen for 
light emitted by a spheroidal bulge. On the near part of the 
galaxy we assume that all the light is extinguished (i.e. the 
bulge is all behind the screen). On the far side the situation 
is reversed, with the bulge un-extinguished. These assump- 
tions are justified by the high inclination of NGC7331 and 
by the luminous bulge, that has an effective radius similar to 
the ring (31 arcsec; Boroson 1981). Within this model, B-K 
on the far side can be regarded as the intrinsic colour of the 
bulge, and the difference between the two half as the colour 
excess due to dust extinction. From the difference in B-K of 
0.5 mag we can thus derive a optical depth in the V-band 
of ~ 0.4, if a galactic extinction law is assumed. In reality 
the dust distribution would have a finite thickness; besides 
scattering and dumpiness should be considered. Inclusion 
of these features would go in the direct ion of reducing the 



extinction for a given quan tity of dust (Boisse & Thoraval 



for the :olour asymmetry of these outer regions. As a check, 
we repeated the correlation using pixels of 15 arcsec, but the 
conclusions are the same: our results are thus independent 



199^ ; [Witt fc Gordon 1996fl . 

Submm images permit us to derive the V-band optical 
depth of the dust distribution in a way independent of ge- 
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ometry. The optical depth tv can be related to the flux at 
A =850-pm using the formula (Hildebrand 1983) 



TV = 



Bx(T) 



where T is the temperature of the dust, I\ is the inten- 
sity (i.e. the flux per beam divided by the integrated beam 
size), B\(T) is the Planck function and the emissivity in- 
dex. Quv is the extinction efficiency in the ultraviolet range 
(1500-3000A), while Q\ is the emission efficiency at a refer- 
ence wavelength Ao- The ra tio Quv /Q\ is q uite u ncertain: 
but from Hildebrand (1983) and from Casey (1991) we have 



derived a mean value Quv/Q\ a ~ 1200 at Ao = 125/xm for 
— 2. The value of Quv/Q xn for = 1.5 has been roughly 



derived from data in Casey (1991): we use Quv/Q\ = 2100 



at Ao = 125 fim. Typically, values of Quv/Q\ Q have un- 
certainties of a factor o f 2-3. The f actor 1/2 comes from 
assuming tuv/tv = 2 (Casey 1991). For the ring (mean 



flux of 70 Jy beam ) we obtain face-on optical depths 
of tv = 2.4 and 9.3, for = 1.5 and 2, respectively. At 
the la isophote tv = 0.4 and 0.9. Uncertainties in tv are 
quite large, because of both the uncertainties in T and in 
the value Quv/Qx - Using temperatures derived with the 
help of KAO data, optical depths are tv — 2.6 a nd 5.4, for 
— 1.5 and 2, respectively. Telesco et al. (1982) derived a 
tv ~ 3 for the reddest part of the galaxy. 

The mass column density of dust can be derived from 
the optical depth. Adopting a mean dust grain radius of 
O.T/im, a density of 3 g cm' 3 and an extinction effic iency 
in the V-band of 1.5 (jffildebrand 1983|; |Casey lggjh, we 



N d = 1.1 M pc" 2 for = 1.5 and 2., respectively (N d = 
0.04 M pc" 2 and N d = 0.1 M pc" 2 at la level). From 
CO ma ps of the east and west side of the ring, Tosaki & 
Shioya (1997) have deduced a mean mass column density 
for the molecular gas of 150 M@pc~ 2 (we note here that 
there are possible uncertainties related to the smaller beam 
size of their observations). An atomic gas co lumn density 



of 20 Mqpc has been derived from Bosma (1981) for the 



position of the ring. The total gas column density for the 
ring is thus N g = 17OM0pc -2 . The gas-to-dust mass ratio 
is then 150 for = 2 (570 for = 1.5). The value derived 
for /3 = 2 is much closer to the galacti c value than IRAS 
based ratios ( Devereux fc Young 199 ^ ) and in ag reement 
with ISO based determination ( [Alton et al. 1998a ; Davies 
et al. 1 39S ) , as well as with optical-NIR radiative transfer 



models (Xilouris et al. 1997, 1998). 

The dust mass inside 150 arcsec is Md = 4 x 10 7 Mq 
for = 2 (M d =1 x 10 7 Mq for = 1.5). From Alton et 
al. (1998a) data, covering a larger area (w 10 arcmin), we 
can derive a dust mass of Md = 1.2 x 10 s Mq (for both 
values of 0). Therefore the total mass is at least three times 
larger than the one we have derived here. This shows that 
there is a large and diffuse component of dust (> 70 per 
cent in mass) that is not associated with the molecular ring. 
It is interesting to note that an extended dust distribution 
may be the reason why some millimetre and submm obser- 
vations J Eales, Wynn- Williams fc Duncan 198S| ; Clement: 



et al. 1393) have failed to detect cold dust, because of the 
small size of of the bolometers us ed compared to the extent 
of the dust ( |Clements et al. 1993 ) . 
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